Stable quark stars beyond neutron stars: can they 
account for the missing matter ? 



■ Peter Minkowski and Sonja Kabana| 

■ f Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, 3012 Bern, 
Switzerland 

E-mail: mink@itp.unibe.ch 



o 

(N 



O 



5— I ' 

O | Laboratory for High Energy Physics, University of Bern, Sidlerstrasse 5, 3012 Bern, 

| Switzerland 

■ E-mail: sonja. kabana@cern.ch 

Abstract. The structure of a spherically symmetric stable dark 'star' is discussed, 
' at zero temperature, containing 1) a core of quarks in the deconfined phase and 

f^) . antileptons 2) a shell of hadrons in particular n, p, A and £~ and leptons or antileptons 

and 3) a shell of hydrogen in the superfluid phase. If the superfluid hydrogen phase 
goes over into the electromagnetic plasma phase at densities well below one atom / 
, (10 fm ) 3 , as is usually assumed, the hydrogen shell is insignificant for the mass and 

the radius of the 'star'. These quantities are then determined approximatively : mass 

■ = 1.8 solar masses and radius = 9.2 km. On the contrary if densities of the order of 
one atom / ( 10 fm ) 3 do form a stable hydrogen superfluid phase, we find a large 

pi 1 range of possible masses from 1.8 to 375 solar masses. The radii vary accordingly from 

9 to 1200 km. 

>• , 1. Introduction 

><■ 

Many authors have found solutions for strange quark stars with shells of hadronic 
composition [1]. We investigate here the structure of a spherically symmetric stable 
dark 'star' as outlined in the abstract. We use for the quark core the Tolman solution 
VI [2] with infinite central energy density and pressure. We also introduce a superfluid 
hydrogen shell which can be very large depending on its inner energy density. For normal 
hydrogen densities the hydrogen shell is very thin and negligible in mass. Objects with a 
mass of few solar masses can in principle be detected by gravitational lensing. However, 
in case of abnormally large initial energy density of the hydrogen shell, we arrive at very 
large masses of the star, which can not be easily detected by lensing. In this case, such 
stars are candidates for 'nonbaryonic' dark matter. The three shells are discussed in the 
following section. For details of this analysis see reference [3]. 

2. The shells of the star 

We consider 3 shells of radii rO (quark shell), rl (hadron shell), and r2 (hydrogen shell). 
The quark core of the star is described according to Tolman's solution VI [2] and energy 
density and pressure are given by 

n - 3 1 j _ 1 1 (1-9 (B/A) r) 

l Je ~ 56 7T GnT 1 dI1U y ~ 56 n GnT 1 (1-(B/A) r) ' 
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Figure 1. 

Temperature in GeV as a function of the baryochemical potential in GeV, 
showing the path followed by the early universe after the QCD phase transition. 

Where A, B are integration constants. We consider the core made up of u,d,s quarks 
and antileptons. allowing for the three chemical potentials associated with charge (hq), 
baryon number (/i#) and total lepton number (fJ,i). The hadronic phase is composed of 
baryons, leptons and antileptons. We use in the hadronic phase the vacuum pressure 
from [4]. The TOV equations allow for a cosmological term, which in the hadronic phases 
must be chosen to cancel the positive vacuum pressure. It thus appears with negative 
sign in the deconfmed core phase. The Gibbs equal pressure condition cannot be met, 
at zero temperature, neglecting a repulsive and very large second virial coefficient in the 
hadronic phase. Thus the equilibrium condition corresponds to equal energy density [3]. 
The hydrodynamic and QCD equilibrium conditions between the two phases (hadronic 
and plasma) together with the condition of locally vanishing overall electric charge 
determine the core radius, mass and all chemical potentials for charge, baryon number 
and lepton number. We derive the mass M and the radius rO of the core: M (quark core) 
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Figure 2. 

Relative fluctuations of the net baryon number normalized to a volume of 1 fm 3 
as a function of the baryochemical potential, following the path of figure 1 

= 0.457 solar masses and r0(quark core) = 3.149 km. The baryon number of the core 
is N b core ~ 0.287 N B sun . Therefore the core itself is not bound, but stabilized 
by the outer hadron shell. We simplify the hadron shell composition structure by 
considering only neutrons, protons and electrons, though near the QCD transition zone 
also E~, A, e + and all antineutrinos intervene. We obtain : M ( hadronic shell ) = 
1.344 solar masses and rl ( hadronic shell ) = 9.047 km. Therefore the total mass of 
both phases inside the radius rl is: M ( total inside rl ) = 1.801 solar masses and the 
baryon number inside r=rl is 1.853 NB(sun). The nonrelativistic approximation of the 
TOV equation used here overestimates somewhat both mass and radius. We proceed to 
calculate the mass and radius of the hydrogen shell which starts at radius rl and ends 
at radius r2. We use here the original approximate relations for superfluid hydrogen 
|, or alike boson-condensed non relativistic system, first described by N. N. Bogolubov 

t We do not include any further light elements, fermions or bosons, in the hydrogen shell. 
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[5]. The results dependent on the parameter H, which is the initial energy density 
of hydrogen at r=rl, normalized to am e m p m 2 H . The total radii, masses and baryon 
numbers are shown in table 1 together with H. The normal hydrogen shell corresponds 
to the first line. 



H 


Radius of dark 'star' 
(km) 


Mass of dark 'star' 
(solar masses) 


N B 
(N B (sun)) 


5 lir 3 


9.204 


1.801 


1.853 


0.290 


411 


2.191 


2.243 


0.295 


674 


3.179 


3.232 


0.310 


1061 


12.63 


12.68 


0.50 


1203 


168.13 


168.18 


0.625 


1208 


271.06 


271.11 


0.75 


1210 


374.02 


374.07 



Table 1. Mass and radius for the complete 'star' with a quark (anti)lcpton core, a 
hadron shell and a hydrogen shell. H is the energy density of the hydrogen shell at 
r=rl normalised to p el — otm e m p m 2 H . 

3. The big annihilation 

We follow the path of the early universe right after the QCD phase transition which 
is important for the genesis of quark stars. In particular fig. 1 shows the temperature 
as a function of the baryochemical potential. At about T ~ 30 MeV the annihilation 
of matter and antimatter begins. We denote this 'the big annihilation'. Note that 
laboratory experiments measured a copious production of antibaryons at a temperature 
of ~ 170 MeV see e.g. [7]. The statistical fluctuations of net baryon number normalized 
to a volume of 1 fm 3 following this path, are shown in fig. 2. They show maxima at 
the phase boundary and at temperatures below 10 MeV. The rest of the condensation 
of such a star needs large gravitational fluctuations too. 

4. Conclusions 

We investigated the structure of a spherically symmetric stable dark 'star'. For a 
low density hydrogen shell we find for the mass and the radius of the 'star' mass = 
1.8 solar masses and radius = 9.2 km. For a high density hydrogen shell supporting 
densities of up to one atom / ( 10 fm ) 3 we find a large range of possible masses from 
1.8 to 375 solar masses. The radii vary accordingly from 9 to 1200 km. Such objects 
are candidates for dark matter of any kind, despite their clearly baryonic nature, if 
they have been formed in such a way that they do not affect nucleosynthesis and the 
cosmic microwave background fluctuations at decoupling. We have calculated the path 
of the early universe after the QCD phase transition in the (T, //#) plane, and the net 
baryon number fluctuations in the absence of gravitational fluctuations. We found that 
the statistical net baryon number fluctuations are large at the QCD phase transition at 
Hb ~ and at T ~ 10 MeV and less and at high \x B . These are regions where such 



5 



stars could be formed. The mass dominating part of the hadron shell is very similar to 
a neutron star. This can be seen in our estimate of the mass : ~ 1.34 solar masses, 
slightly below the Chandrasekhar bound. The distinctive feature to a neutron star is 
that the hadron shell does not extend to the center of the 'star'. The quark core taken 
for itself and according to Tolmans universal solution, with p e oc r~ 2 towards the 
center is not stable. This means that it contains less nucleons counted by number than 
counted by mass. Yet this does not imply any instability of the 'star' as a whole, as it 
would be the case for a one-phase neutron star. The phase boundary between quark 
core and hadron shell with its surface pressure stabilizes the 'star'. Despite genuine 
baryonic origin, besides the (anti)leptonic parts, such 'stars' are indeed candidates for 
any kind of dark matter. If 1.8 solar masses indeed characterize the mass of these 'stars' 
then they could be observed through gravitational microlensing of the galactic halo 
towards the Magellanic clouds [6]. However masses much larger than few solar masses 
are hardly detectable by lensing. We also note, that the composition of the 'star' 
considered here lacks heavy elements and hence the metallicity nearby is unusually low 
compared to burning stars of similar mass. They could possibly emit and absorb H lines, 
e.g. Lyman — a lines taking into account gravitational redshift (which is z=0.56 for a 
mass of 1.8 solar masses) and high pressure compared to normal conditions. Another 
possibility is that they are part of a gravitating system including one or more shining 
stars. 
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